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Abstract 



Reduction is a process that uses symmetry to lower the order of a Hamiltonian 
system. The new variables in the reduced picture are often not canonical: 
there are no clear variables representing positions and momenta, and the 
Poisson bracket obtained is not of the canonical type. Specifically, we give two 
examples that give rise to brackets of the noncanonical Lie-Poisson form: the 
rigid body and the two-dimensional ideal fluid. From these simple cases, we 
then use the semidirect product extension of algebras to describe more complex 
physical systems. The Casimir invariants in these systems are examined, 
and some are shown to be linked to the recovery of information about the 
configuration of the system. We discuss a case in which the extension is not a 
semidirect product, namely compressible reduced MHD, and find for this case 
that the Casimir invariants lend partial information about the configuration 
of the system. 
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This paper explores the Casimir invariants of Lie-Poisson brackets, which generate the 
dynamics of some discrete and continuous Hamiltonian systems. Lie-Poisson brackets are a 
type of noncanonical Poisson bracket and are ubiquitous in the reduction of canonical Hamil- 
tonian systems with symmetry. Casimir invariants are constants of motion for all Hamil- 
tonians; they are associated with the degeneracy of noncanonical Poisson brackets. Finite- 
dimensional examples of systems described by Lie-Poisson brackets include the heavy top 
and the moment reduction of the Kida vortex, while infinite-dimensional examples include 
the 2-D ideal fluid, reduced magnetohydrodynamics (MHD), and the 1-D Vlasov equation. 
(See Ref. |T| and references therein for a full review.) The Casimir invariants determine the 
manifold on which the system is kinematically constrained to evolve. Understanding the 
nature of these constraints is thus of paramount importance. 

In Section | we examine specific Lie-Poisson brackets, namely, those that arise from the 
reduction to Eulerian variables of a Lagrangian system with relabeling symmetry. We make 
use of two prototypical examples, the rigid body (finite-dimensional) and the 2-D ideal fluid 
(infinite-dimensional), and we interpret their Casimir invariants. This is done to motivate 
the introduction of such brackets and to show their physical relevance. In Section |I| we turn 
to building Lie-Poisson brackets directly from Lie algebras by the procedure of extension. 
We introduce the semidirect product extension and illustrate it with two physical examples: 
the heavy top and low-beta reduced MHD. Finally in Section [TTT] we look at a nonsemidirect 
example, compressible reduced MHD, and discuss work in progress. 



I. LIE POISSON BRACKETS AND REDUCTION 

For our purposes, a reduction is a mapping of the dynamical variables of a system to a 
smaller set of variables, such that the transformed Hamiltonian and bracket depend only on 
the smaller set of variables. (See for example Ref. [|for a detailed treatment.) The simplest 
example of a reduction is the case in which a cyclic variable is eliminated, but more generally 
a reduction exists as a consequence of an underlying symmetry of the system. We present 
two examples of reduction. 



A. Reduction of the Free Rigid Body 

The Hamiltonian for the free rigid body is an unwieldy function of three Euler angles 
<f>,ip,9 and their conjugate momenta p<f>,p^,P8. The motion is described by Hamilton's 
equations using the canonical bracket 

_ df dg dg df | df dg dg df | df dg dg df 

can d<p dp^ dcj) dprf, dip dpip dip dp^ dd dpe dd dpg 

Here we have 3 degrees of freedom (6 coordinates), the configuration space is the rotation 
group SO (3), and the phase space is its contangent bundle T*SO(3). 

A reduction is possible for this system. In terms of angular momenta £{ about the 
principal axes, we have 

i=i * 
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Under this (noncanonical) mapping, the bracket obtains the Lie-Poisson form 



if..}- 

The equations of motion generated by this bracket and H(£i, £ 2 , are Euler's equations 
for the rigid body, 

&i j j &j ^k i 

I j lk 

where i, j, k are cyclic permutations of 1, 2, 3. The energy is conserved, and so is the quantity 
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which commutes with any function of i. Such a functions are called Casimir invariants (or 
Casimirs for short). Casimirs are conserved quantities for any Hamiltonian, so they tell us 
about the topology of the manifold on which the motion takes place. For the simple case 
of the rigid body, the motion takes place on the two-sphere, S 2 (not in physical space, but 
in angular momentum space). The symmetry that permits the reduction is the invariance 
of the equations of motion for (0, 6 l ,^,p^,pe) under rotations (elements of 5*0(3)). This 
symmetry amounts to the freedom of choosing axes from which the Euler angles are mea- 
sured. In that sense it is a relabeling symmetry, since the choice of axes amounts to making 
"marks," or labels, on the rigid body. 

We shall say that the original system is a Lagrangian description (by analogy with 
the fluid case below) because at any time the exact configuration of the system (including 
orientation) is known, whereas the reduced system is Eulerian because only the angular 
momentum of the body is known. 



B. Reduction of the 2 D Ideal Fluid 

As our prototype for reduction in an infinite dimensional system we take an ideal 2- 
D fluid confined to some domain, D. The Lagrangian description involves fluid elements 
labeled by some coordinate a, which is usually taken to be the initial position of the fluid 
elements. These labels are analogous to the choice of axes in the rigid body example above 
(the "marks" on the rigid body). The Hamiltonian functional is 



H[q-M 



D 



7T 



O P(M) 



dq 
da 



d a, 



where q(a,t) is the position of the fluid element labeled by a and ir(a,t) is its momentum. 
The Jacobian of the transformation from the labels a to the position of the fluid elements at a 
later time is \dq/da\. The density p is taken to be constant, and the pressure p(a, t) appears 
here as a Lagrange multiplier that enforces the incompressibility condition, \dq/da\ = 1 (see 
Ref. H). This Hamiltonian together with the canonical bracket 
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'5F5G SG SF' 



d a 



generates the equations of motion for a fluid in Lagrangian variables. The information about 
the position of every fluid element at any time is contained in the model. The dynamical 
evolution of the system is independent of the particular choice of labels for the fluid elements. 
This relabeling symmetry of the initial condition labels, a, suggests a reduction. 

We introduce the streamfunction defined by v (x, t) = (—d y 0, d x 0), so that V • v — is 
automatically satisfied, and the vorticity u;(x, t) = V 2 0- The mapping from the Lagrangian 
momentum to the Eulerian velocity field is 



Pov(x,t) 



D 



7r(a, t) <5(x — q(a, t)) d 2 a . 



We take p — 1 for simplicity. Taking the curl of v and dotting with z gives 

c<j(x, t) = z ■ V x t>(x, t) — J €ij 7Tj(a, t) dj5(x — q(a, t)) d 2 a , 

where repeated indices are summed, d{ := d/dxi, and the antisymmetric symbol is defined 
by €12 = 1. The variation of is 

5^(x,t) = (e ij 57r i (a,t)d j 5(x- q(a,t)) - e^ir^t) djd k 5(-x- q(a,t)) Sq k (a,t)) d 2 a, 



which we can insert into 

SF 
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x=q(a,t) 



SF r SF 
"l - = ~ / e ij ^l ' *) <9jV5(x - q(a, t)) — d 2 x. 
Sq JD Suj 

We then insert these two expressions into the canonical bracket. After some manipulation 
involving integration by parts (we assume boundary terms vanish) we obtain the Lie-Poisson 
bracket 



{F,G}= f 

JD 



SF SG 

Suj ' Suj 



d 2 x , 



where 



[/,<?] 



df_ dg_ _ d£ dg_ 

dx dy dy dx 



Note that the incompressibility of the fluid (the fact that the Jacobian \dq/da\ is unity) was 
not used in the derivation of this bracket. However, in order to write the Hamiltonian in 
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terms of u one must introduce the streamfunction 0, which is possible only if V • v = 0. The 
equation of motion generated by the bracket and the transformed Hamiltonian 

H[u] = --f <puod 2 x = \f |V0| 2 d 2 x 

is just Euler's equation for an the ideal fluid 

w(x) = -[0,w] . 

This has a Casimir given by 

C[u]= I f(u(it))d\ 
Jd 

where / is an arbitrary function. The interpretation of this invariant is given in detail in 
Ref. |j. It implies the preservation of contours of a;, so that the value uo on a contour labels 
that contour for all times. This is a consequence of the dissipationless and divergence-free 
nature of the system. Substituting f(uS) = u n we also see that all the moments of vorticity 
are conserved. By choosing f(uS) = 6 (no (it) — u ), a heavyside function, it follows that the 
area inside of any cu-contour is conserved. 

II. EXTENSIONS AND THE SEMIDIRECT PRODUCT 

We now investigate systems involving Lie algebras by extension, a procedure for com- 
bining two or more Lie algebras to make a new Lie algebra. There are a myriad of ways to 
extend algebras, and we will only touch on a few here. All the extensions discussed here 
have their equivalent for Lie groups, but we choose the algebra approach here because it 
leads more directly to a Lie-Poisson bracket. In this section we let f) be an extension of the 
Lie algebra g by the algebra t>. The elements of f) are written as 2-tuples, (£, rj), where £ G g 
and rj e d. 

The simplest extension is the direct product of Lie algebras. Let £ and £' be elements of 
a Lie algebra g and rj and rf be elements of a vector space o (which is an Abelian Lie algebra 
under addition). The direct product of these two algebras is an algebra f) with bracket 

K^),(£W) ]:=([£,£'], [WD • 

Given the same g and d as above there is a less trivial way to make a new Lie algebra 
called the semidirect product with an operation defined by 

i(t,v),(e,v')]--=([t;,e]ALv'] + [v,e])- 

A simple example of a semidirect product structure is when g is the Lie algebra so(3) 
associated with the rotation group 5*0(3) and o is M 3 . Their semidirect product is the 
algebra of the 6-parameter Galilean group of rotations and translations. Both the elements 
of g and o can then be represented by vectors in M 3 , with bracket [£,77] = £ x 77, the cross 
product of vectors. Since f) is itself a Lie algebra, it can be extended again as needed to 
make an n-fold extension (an algebra of n-tuples). 
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We can build Lie-Poisson brackets from these algebras by extension.Q For an n-fold 
extension i) of the Lie algebra g, we define 

where \x G i)*, the dual of f) under the pairing (,):Fj*xf)— >R. The dynamical variables of the 
system are the elements of the n-tuple \x = fi{t). These elements may be fields or variables, 
so the Lie-Poisson bracket derived from an algebra by extension generates the dynamics for 
a system involving several dynamical quantities. The functions (or functionals) F and G 
are maps from fj* to R. Here 5/5fi is a derivative or a functional derivative, depending on 
the dimensionality of the algebra (finite or infinite). For n = 1, Fj = so(3), we have // := £ 
and we recover the bracket for the free rigid body (Section [I A|) . The overall sign of the 
Lie-Poisson bracket has to do with left- or right-invariance of vector fields and will not be 
discussed here (See Ref. 0). 

Using this procedure to make a Lie-Poisson bracket from a direct product of algebras 
leads to a sum of n independent brackets. This will not interest us further, since the 
coupling between dynamical variables can only come from the Hamiltonian. However there 
are interesting physical examples of a direct product structure. (This is the case for the 
model in Ref. |5], although a coordinate transformation is needed to exhibit the structure.) 

We illustrate the process of building a Lie-Poisson bracket from a semidirect product of 
algebras by two examples, which are extensions of the rigid body and ideal fluid examples 
of Section |I[ 



A. The Heavy Top 

The Lie-Poisson bracket for the semidirect product of the rotation group SO (3) and the 
vector space R 3 is 

ri , „ (df dg\ fdf dg Of dg\ 
where a denotes a 3-vector. By using 

3 p2 

2=1 1 

where c is a vector representing the position of the center-of-mass, we get the prototypical 
example of a semidirect product system, the heavy rigid body (in the body frame): 

• i,- - Ik „ „ .4 atj in a k 



Ej 4 + ctj c k - OL k cj , cxi 



Ij Ik h I 



J 



where i,j,k are cyclic permutations of 1,2,3. The vector a rotates rigidly with the body, 
which is always true for a Hamiltonian quadratic in £. The Casimirs for this bracket are 

d = a 2 , C 2 = I ■ a . 
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Looking at the bracket as derived by reduction of the heavy top in Euler angles (as we did 
in Section IA, but here with gravity), the Casimir C 2 expresses conservation of p^, since 
is cyclic. Knowing a does not lead to a determination of the orientation of the rigid body: 
there is still a symmetry of rotation about a. Taking the semidirect product has led to the 
recovery of some of the Lagrangian (configuration) information. 



B. Low-beta Reduced MHD 



The semidirect product bracket for two fields uj and i/j is 



{F,G}= f 

JD 





f 

UJ 

i. 


~5F 


5G~ 




~5F 5G~ 


+ 


~5F 


5G~ 




L 


5u) 


5u) 


5uj ' 5tp 




5uj 





d 2 x 



If uj = V 2 0, where <ft is the electric potential, ip is the magnetic flux, and J = V 2 ip is the 
current, then the Hamiltonian 



HM\= l -J D (|V0| 2 + |V^| 2 ) d' 
with the above bracket gives us 



x 



UJ 



[u, 0] + [ip,J] , V = 0] > 

a model for low-beta reduced MHD.I (Ref. |] contains a system with a similar structure, but 
for waves in a density-stratified fluid.) 

The bracket has two Casimir invariants, 



D 



f(iP)d 2 x, C 2 [u;;V] 



D 



uj g(ijj) d 2 x. 



The first has the form of the Casimir for 2-D Euler of Section 
pretation. To understand the second one we let g(ip) = 8(ip — ipo 
this case we have 



and has the same inter- 
a heavyside function. In 



C 2 [uj-ij] 



*0 



uj d 2 x, 



where \l/o represents the (not necessarily connected) region of D enclosed by the contour 
i/j = -00, and d^o is its boundary. The contour d^o moves with the fluid, so this just 
expresses Kelvin's circulation theorem: the circulation around a closed material loop is 
conserved. 

This theorem is true for any -0-contour, therefore it holds in the region between two 
contours ip — ipo an d ip = ipo + 5. Letting 5 — > we see that the two contours delineate a 
"line" of fluid elements with value tpo of the magnetic flux. Knowledge of the value of ijj on a 
fluid element thus only determines which contour it is on, but not its location on the contour. 
Therefore, there is still a relabeling symmetry: the fluid elements can be shifted around the 
contour without changing the Casimirs C\ and As with the heavy top, the semidirect 
product has led to the recovery of some, but not all, of the Lagrangian information. 
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C. Putting Labels on a Rigid-body 



Remember that taking a semidirect product restricted the symmetry group of the body 
to rotations about a. If we take another semidirect product to get 

r* n „ fdf dg\ fdf dg df dg\ a fdf dg df dg\ 

where (3 is a 3-vector, we have a bracket that can model a rigid body with two forces acting 
on it, for example a charged, rigid insulator in an electric field. The new bracket has Casimirs 

d = a 2 , C 2 = (3 2 , C 3 = a-(3. 

The angular momentum i has disappeared from the Casimirs. This is because knowing a 
and f3 completely specifies the orientation of the rigid body (unless the two are colinear). In 
other words, by taking semidirect products we have reintroduced the Lagrangian information 
into the bracket. Note that taking more than two semidirect products is redundant as far as 
the Lagrangian information is concerned: knowing the orientation of more than two vectors 
does not add new information. This is reflected by the fact that the number of variables 
minus the number of Casimirs is six for two or more "advected" quantities. This is the 
dimension of T*SO(3), the original phase space (before reduction). 



D. Advection in an Ideal Fluid 

We now take a second semidirect product for the ideal fluid, say low-beta MHD with a 
second advected quantity, the pressure p. In that case we get a model for high-beta reduced 
MHD.i The Casimir is 

C[ip;p] = / f(ijj,p)d 2 x, /arbitrary. 

JD 

This Casimir amounts to being able to label two contours. Locally this permits a unique 
labeling of the fluid elements as long as p and ip are not constant in some region. However, 
globally there is some ambiguity, because contours can cross in several places. Thus, in 
the infinite-dimensional case the semidirect product is not equivalent to recovering the full 
Lagrangian information, unless the contours do not close, are monotonic, and nonparallel 
(Vip x Vp does not vanish). A third advected quantity will in general break this degeneracy. 
Note that if the advected quantities label the fluid elements unambiguously at t — then 
they will do so for all times. 



III. BEYOND THE SEMIDIRECT PRODUCT: COCYCLES 



In general there are other ways to extend Lie algebras besides the semidirect product. 



One example is the model derived in Refs. [10] and |TT] for 2-D compressible reduced MHD. 
The model has four fields, and is obtained from an expansion in the inverse aspect ratio of 
the tokamak. The Hamiltonian is 



S 



H[u,v,p,< 




d 2 x, 



where v is the parallel velocity, p is the pressure, and (3 is a parameter that measures 
compressibility. The bracket is 



{A,B} 




The term proportional to (3 is an obstruction to the semidirect product structure, and it 
cannot be removed by a coordinate transformation. 

The theory that deals with the classification of extensions is Lie algebra cohomology. In 
general the way to extend a bracket is by adding a nontrivial cocycle. Though a priori there 
are an infinite number of ways to make an extension, for low dimensions, after allowing for 
coordinate transformations, very few possibilities remain; we have classified these in Ref. [12. 
We have also found all the Casimir invariants for the low-dimensional brackets (five fields 
or less). 

The Casimirs of the above bracket are 



CM = 



D 



fWd\ 
v h(ip) d 2 x, 



D 



D 



C 2 [p;^] = l_ pg(*P)d x 
Ci[uj,v,p,ip) 



D 



d x. 



Finding the invariant C4 directly from the equations of motion would be tedious, but is 
straightforward from the bracket. These Casimirs do not allow a labeling of the fluid ele- 
ments. The meaning of invariants of the form of Ci, C2, and C3 was discussed in Sections [T~B 



and |11B| : the total magnetic flux, pressure, and parallel velocity inside of any ^-contour are 
preserved. To understand C4 we use the fact that uj = V 2 </> and then integrate by parts to 
obtain 



C 4 [uJ,v,p, ip] 



■ W + y ^ k'(ip) d 2 x. 



The quantity in parentheses is thus invariant inside of any ^-contour. It can be shown that 
this is a remnant of the conservation in the full MHD model of the cross helicity, 



V 



v ■ B d x . 



D 



at second order in the inverse aspect ratio, while C3 is a consequence of preservation of this 
quantity at first order. Here B is the magnetic field. As for C\ and C2 they are, respectively, 
the first and second order remnants of the preservation of helicity, 

W = [ A B d 2 x, 
Jd 

where A is the magnetic vector potential. 
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IV. CONCLUSIONS 



We gave an introduction to the reduction of physical systems based on their symmetries. 
The prototypical examples were shown, the rigid body and the 2-D ideal fluid. For these 
two cases some information about the configuration of the system was lost after reduction, 
correponding to the symmetry used to reduce the system. 

The semidirect product allowed us to build larger brackets from a "base" algebra in a 
systematic manner. We were thus able to describe the heavy top and low-beta reduced 
MHD. Examining the invariants, we concluded that the semidirect product had recovered 
some or all of the Lagrangian information. 

For general extensions (not necessarily semidirect) things are different: the Lagrangian 
information is not necessarily a consequence of the Casimirs. However, for compressible 
reduced MHD the Casimirs represent constraints that are remnants of invariants of the full 
MHD equations from which the model is derived asymptotically. 

As mentioned in Section |TTT| , when considering a general extension, all brackets can be 
reduced to a small number of normal forms, at least for low-dimensional extensions. It will 
be interesting to see if physical systems can be found that are realized by these brackets, both 
in the finite and infinite degree-of-freedom cases. We are currently investigating a toy model 
that we call the Leibniz top, which is one of the simplest non-semidirect system one can 
build. It is a straightforward generalization of the Lagrange top (a heavy top with Ii = I 2 ). 
The Lagrange top is integrable, and we have found that so is the Leibniz top. 
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